arXiv:l506.0663lv3 [hep-ph] 3 Nov 2015 


Extending theories on muon-specific interactions 

Carl E. Carlson 1 : E| and Michael Freid 1,2 ’[0 

1 Nuclear and Particle Theory Group, College of William and Mary, Williamsburg, Virginia 23187, U.S.A. 

2 Thomas Jefferson National Accelerator Facility, Newport News, Virginia 23606 

(Dated: August 7, 2015) 

The proton radius puzzle, the discrepancy between the proton radius measured in muonic hydro¬ 
gen and electronic hydrogen, has yet to be resolved. There are suggestions that beyond the standard 
model (BSM) physics could resolve both this puzzle and the muon anomalous magnetic moment 
discrepancy. Karshenboim et al. point out that simple, nonrenormalizable, models in this direction 
involving new vector bosons have serious problems when confronting high energy data. The prime 
example is radiative corrections to W —> gv decay which exceed experimental bounds. We show 
how embedding the model in a larger and arguably renormalizable theory restores gauge invariance 
of the vector particle interactions and controls the high energy behavior of decay and scattering 
amplitudes. Thus BSM explanations of the proton radius puzzle can still be viable. 


I. INTRODUCTION 

The 7a discrepancy between the proton radius mea¬ 
sured via muonic hydrogen Lamb shift DO [21, and mea¬ 
sured via electron scattering or electronic hydrogen 
atomic level splittings (including the Lamb shift in or¬ 
dinary hydrogen) [3], remains unresolved. Assuming the 
discrepancy is not due to experimental error, we can con¬ 
sider the explanation that there exists unaccounted for 
physics which masks itself as a difference in proton size. 
Such an effect must lower the muonic hydrogen 2 S state 
by 310 geV to match experiment. There have been sev¬ 
eral proposed beyond-the-standard-model (BSM) theo¬ 
ries that could explain such an effect. See mu for re¬ 
views. 

Furthermore, there exists a discrepancy of 3 or more a 
between the measured muonic g — 2 and that predicted 
by the standard model m- A tantalizing hope and goal 
of these BSM theories is to solve both the proton radius 
and the muonic g — 2 problems simultaneously. 

One idea in the BSM regime is the existence of a new 
force that couples preferentially to muons mm- How¬ 
ever, a number of constraints must be addressed ma¬ 
ns]. Many of them can be avoided if we suppose the 
new force carrier (s) couple only to muons on the leptonic 
side (or at least much more strongly to muons than to 
other leptons), and only to first generation particles on 
the hadronic side, with the hadronic coupling propor¬ 
tional to the electric charge. A neat scheme of this sort 
has been proposed by Batell et al. SB and further stud¬ 
ied by Karshenboim et al. [15| . Our own study will be 
more generic, along the lines of models studied by Rislow 
and one of the present authors [15]; these models in par¬ 
ticular have smaller electron couplings than in mm 
and a different pattern of parity violation such that cor¬ 
rections to the muonium hyperfine splitting and 133 Cs 
weak charge will not arise. 
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The new muon coupling will affect other processes in¬ 
volving muons, and one immediately thinks of the muon’s 
anomalous magnetic moment, (g — 2)^, with its circa 
2.1 ±0.7 part-per-million (ppm) discrepancy between the 
standard model and experiment. In the present context, 
it was early noted [TO] that if the mass of the BSM force 
carrier is very light, its effect on (g — 2) M can explain the 
discrepancy directly. If the mass of the new force carrier 
is not very light, one has to arrange a fine-tuning involv¬ 
ing at least one more BSM force carrier to keep the effect 
of the new physics on (g — 2)^ at the correct size. This is 
by now well understood mm, and one can even turn 
it into a positive by solving the (g — 2)^ problem and the 
proton radius puzzle with the same model. 

Further, if the BSM force carrier is not too heavy, it 
can contribute as a bremsstrahlung or radiative correc¬ 
tion to decays involving muons. For example, the decay 
K —> gv(f), where <f> is the new force carrier was studied by 
Barger et al. m for the case where <j> does not visibly de¬ 
cay. If there is a small coupling allowing </> —► e + e _ , then 
the e + e _ spectrum would show a visible bump above the 
QED background in K — >• gve + e~ , as was considered 
in [151 . 

Karshenboim et al. na called special attention to the 
need for embedding any new vector force carrier into a 
renormalizable theory, or at least into a theory where 
the scattering or decay amplitudes do not grow with en¬ 
ergy when the energy is high. They did this by studying 
radiative corrections to W —> gu decay. The width of 
this decay is unambiguously predicted in the standard 
model, and is measured with an uncertainty of about 
2%. Karshenboim et al. showed that if there be just a 
low mass vector boson that coupled to muons (but not to 
IT’s or i/’s), with a coupling chosen to explain the energy 
splitting deficit in muonic hydrogen, then its contribu¬ 
tion T(W —> gvV) would exceed the measured width of 
the IF by a large margin. 

However, in a well behaved and renormalizable the¬ 
ory, the growth of amplitudes with energy cannot go 
unchecked. Unitarity imposes limits on the energy be¬ 
havior of scattering amplitudes, and if one is using the 
conventions of (say) Bjorken and Drell [T7] or of Peskin 
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and Schroeder IB] , an amplitude in a single partial wave 
must not grow with energy at high energy ( i.e if the 
amplitude grows like (energy)" for large energies, then 
n < 0). Nonrenormalizable theories are known for their 
ultraviolet divergences in loops, but their excessive en¬ 
ergy dependence can also appear at tree level in the form 
of unitarity violations. A known historical example is 
the amplitude for v e v e —> W + W~ in a simple vector bo¬ 
son theory m ■ The calculation from just diagram []J a) 
gives an amplitude that is asymptotically in a single par¬ 
tial wave that grows like E 2 as the center-of-mass energy 
E —> oo. The Weinberg-Salam extension of the theory 
also has a Z-boson diagram, 0b), which is significantly 
smaller thanjlja) at threshold but asymptotically cancels 
the offending energy behavior and restores perturbative 
unitarity [ 201 . A general study by Llewellyn Smith has 
shown that the need to satisfy unitarity bounds leads to 
a Yang-Mills structure for many theories involving vector 
bosons [2T L 



FIG. 1: The illustrative process vv —> W + W . 

In this paper, we consider new vector and (when 
needed) axial-vector bosonic interactions that couple to 
the muon and the proton but do not couple or couple 
weakly to the electron and most other particles. Again, if 
this is all we have, the result of Karslienboim et al. shows 
that the region of parameter space which solves the pro¬ 
ton radius problem does not occur in the allowed param¬ 
eter space given by the known decay of the W. Inspired 
by (2UJ [21] > we add an additional triple boson vertex in 
the Lagrangian, giving an interaction involving the stan¬ 
dard W-boson, the new vector particle 0y, and a fur¬ 
ther vector boson with the same mass as the W. We 
call this newest boson a “shadow W,” denoted W s with 
m(W s ) = mw■ We also include, when needed, a corre¬ 
sponding axial vector triple boson interaction, involving 
the shadow W s , the ordinary W, and the 0 a- The inclu¬ 
sion of the W s makes the 0 interactions gauge invariant 
or current conserving, arguably fixes the nonrenormaliz- 
ability of the original interaction, and, as we shall show, 
definitively pushes the constraints on the couplings due 
to W decay far away from the coupling strength param¬ 
eter region necessary to solve the proton radius problem. 
Thus it can be a plausible candidate for a BSM solution 
to the proton radius problem. 

We note that a current conserving theory with mas¬ 
sive bosons (<j>v and 0a) and shadow W’ s, gives high 
energy results, e.g. for radiative corrections to W decay, 
very much like a theory with a massive scalar boson 0 S 


plus, when needed, a corresponding pseudoscalar boson 
0 p . We briefly display such a scalar theory, and show 
that decays of the W involving such a scalar and pseu¬ 
doscalar do not restrict the necessary parameter space 
needed for solving the proton radius problem with scalar 
exchanges. 

We should also note that though our theory is well- 
behaved and seems likely to be renormalizable (as argued 
by Llewelyn Smith |21jh it is not yet a full theory em¬ 
bedded into the standard model (SM). Further work will 
be required to show how such a theory can be embed¬ 
ded into the SM. For now, we simply consider our theory 
as a phenomenological application of some BSM physics, 
containing features that a full theory must contain and 
controlling the high energy behavior of scattering and 
decay amplitudes. 

In the following, the bulk of our work concerns the 
new vector or axial vector bosons, and is described in 
Sec. |II} We also include some comments on why the cor¬ 
responding radiative corrections to Z —> decay are 

innocuous. Results for the scalar case are given in a short 
Sec. |III[ and conclusions are offered in Sec. [TV] 


II. VECTOR THEORY 

We start with an interaction Lagrangian similar to [32] 
where 0v interacts with a muon (and proton) via the the 
explicit vector coupling Cy (Cy) and where 0a interacts 
with a muon (and proton) through the axial vector cou¬ 
pling (C^). For brevity of notation, it is understood 
that 0 without a subscript represents either 0y or 0a in 
this section. We also include an additional 3-boson inter¬ 
action |21| term involving the 0, the ordinary W, and a 
third boson, with coupling strength equal to Cy (or G^) 
as is necessary to make the decay W —► \iv0 gauge in¬ 
variant. The third boson is the shadow W, denoted W s , 
which couples to the muon in the same manner as the W 
and has m\y g = mw- 

The new interaction terms in the Lagrangian are, 


Cint =~0x + Cy V’jGVj 

- </>a CaVV7 A 75Vv + C^p7 A 75</7 

- iC»e ijk W' a W 3 p d a W k * + i {C^ terms} 

V’m7 A ( 1 -7bM/ W~ x +h.c., (1) 


2y2 

where in the Cy terms, 


^ o W~ , 

C W+ a , 

w^0 r, (2) 


with V —► A for the terms; is the totally anti¬ 
symmetric Levi-Civita symbol. Note that we could use 
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two shadow W's (one vector and one axial vector), but 
appear to need only one. Further note the different signs 
of the Cy and Yang-Mills terms necessary for gauge 
invariance, and that we have included an interaction of 
the W s with the charge changing muon current. 

If the Cy and Cy have the opposite sign then there 
exists an additional attractive force between the muon 
and the proton through the interaction with the <j)y. This 
additional force will create a difference between the 2 S- 
2 P Lamb shift in muonic hydrogen and hydrogen as us¬ 
ing 


iCygw 




iM = (pi+w) 2 7“(i-7 ) 


/ „ (Pl+P2)„(Pl+P2)u 

- / ypv ^72 

-7 m (1-7 5 ) 


(pi +P2) 2 ^ m 2 


w 


x ^0(k+p 3 r+g^(-p 3 + pi +p 2 ) a 
+ 9 av {-Pi-P2-kf)^(p 2 ) (4) 


AE(2S-2P ) 


ICyCy-l rnl(m r a ) 3 
47 t 2 (to 0 + m r a ) 4 


( 3 ) 


where m r is the reduced mass of the (muonic) hydrogen 
system. The contribution to AE[2S-2P) from the axial 
coupling C^ is very small. 

To account for the energy difference that can be inter¬ 
preted as a proton radius difference, there must be an 
extra 310 /xeV in the 2S-2P Lamb shift of muonic hydro¬ 
gen tun, The parameter Cy necessary to satisfy this 
constraint is plotted as the green band outlined by solid 
lines in Fig. [3] where Cy = \Cy\ = Cy. 

Furthermore, the introduction of new <py and (f>A in¬ 
teractions with the muon will shift the muon anomalous 
magnetic moment. The vector and axial vector couplings 
affect the anomalous moment with opposite signs and can 
be tuned to account for the known discrepancy between 
theory and experiment of muonic g — 2 [12]. If Cy is set 
to satisfy the proton radius problem, then the allowed 
region for C 1 ^ from the muon g — 2 constraint is shown 
by the green band outlined by dashed lines in Fig. [3] 

We now move on to consider a constraint emphasized 
by Karshenboim et al. m, that the branching ratio of 
W — > fjbv<j>y plus W —> gv(j)A must be less than 4 per¬ 
cent (twice the error in the W width as measured by the 
Tevatron). Without the inclusion of a 3-boson interac¬ 
tion, this constraint eliminates the region of the (Cy , 171 ^) 
parameter space required to explain the proton radius 
puzzle. This decay is calculated from the Feynman dia¬ 
grams given in Fig. [2] 




From |T]) we can derive the necessary Feynman rules 
to compute this decay amplitude as 


where k is the W 4-momentum, p 3 is the muon 4- 
momentum, p 2 is the neutrino 4-momentum, and p 3 is 
the (fiy 4-momentum. Here we have focused on the vec¬ 
tor contribution to the W decay, but one can easily show 
that the axial vector contribution is equivalent up to an 
overall minus sign (which is irrelevant to the decay am¬ 
plitude squared). 

Letting the muon and neutrino mass be zero, we find 
(to leading order in m^/mw) 


IV = 


G F m 3 w [( C £) 2 + (C ^) 2 


96-\/2' 


7T 


i 2 m W r 1 m W , ” 1 
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m~ 
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37 

y 



( 5 ) 


Keeping the muon mass would only give multiplicative 
corrections to the coefficients like (1 + 0 (m 2 /m^)). 

This decay width has a strikingly different dependence 
on rricf, compared to [15] , The 1 /m 2 ^ dependence found 
there that came from the longitudinal component of the 4> 
polarization is canceled by the inclusion of the W s prop¬ 
agator in the second diagram of Fig. [2] Thus at lead¬ 
ing order in m^/mw, the mass divergence is logarithmic 
and not inverse polynomial. This logarithmic dependence 
pushes the constraints from W decay far away from the 
desired parameter space of 67) and m^. 

The contribution of C 1 ^ in (pj) can be obtained in terms 
of Cy using the constraint from (g — 2) M [T^- The con¬ 
straint from W decay eliminates the region of (Cy,m^) 
above the top curve, the shaded red area, in Fig. [3] The 
values of Cy below this area are allowed by this con¬ 
straint. 

Another constraint on Cy occurs from transitions be¬ 
tween 3 d and 2 p orbitals in muonic 24 Mg and 28 Si m 
□22 El. At two standard deviations, this constraint is 
plotted as the shaded orange area bordered below by a 
solid black line in Fig. [3] where allowed values of Cy exist 
on and below this line. 

Note that an additional constraint due to muonium 
hyperfine splitting discussed in m is not relevant here 
since </> does not couple to the electron (or the coupling 
can be kept quite small). For similar reasons, we do not 
have a constraint on C 1 ^ from a new parity nonconserv¬ 
ing interaction contributing to the weak charge in 133 Cs, 
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FIG. 3: The parameter space necessary to satisfy experimen¬ 
tal constraints. Solid lines refer to constraints on Cy. Dashed 
lines refer to constraints on C^J. The green band, outlined by 
solid lines, is the constraint on Cy necessary to solve the 
proton radius problem (±2a). The shaded red region is the 
restricted region of Cy due to the constraint that the branch¬ 
ing ratio for W goes to pv(j)y + must be less than 4% 

under the assumption that solves the muonic g — 2 prob¬ 
lem. The shaded orange region is the restricted region on Cy 
due to energy splittings in muonic Mg and Si at 2a. The 
green band, outlined by dashed lines, is the constraint on Ca 
necessary to solve the muonic g — 2 problem (±2<r) under the 
assumption Cy solves the proton radius problem (±2<r). 


significantly opening up the allowed parameter space for 

(C^m^). 

In Fig. [3] we see that there are broad regions of param¬ 
eter space for which we can find values of Cy , (7^, and 
m<t, that simultaneously solve the proton radius puzzle 
and the muonic g — 2 discrepancy while satisfying the 
considered experimental constraints. 

For completeness, we also comment on radiative cor¬ 
rections to Z —>■ p + p~ decay, namely Z —> p~g + (j>v and 
p~ p + 4>A decay as represented in Fig. [ 4 J 



+ 


FIG. 4: Z -A g~g + (j> 



This decay amplitude is 


iM=l 2^ w c y 6a{k)e ^ P3Mpi) 


- 7“ ( - 1 + 2ai uH w - I7 5 ) }“<!» 


where k is the Z 4-momentum, pi is the muon 4- 
momentum, P 2 is the anti-muon 4-momentum, and p$ 
is the <pv 4-momentum. As with the W decay, here we 
only focus on the vector contribution to the Z decay, but 
one can easily show that the axial vector contribution is 
equivalent up to an overall minus sign (which is irrelevant 
to the decay amplitude squared). 

In this case, cancellations between the two diagrams 
ensure the Ward identity is satisfied. Therefore, there is 
no poor behavior at high energies when the <f> is longi¬ 
tudinally polarized. This is seen in the logarithmic de¬ 
pendence of the decay width Q on (which resembles 
that of the W decay Q), 


T z = 


G F m 3 z [((7(0 * 2 + (C^) 2 ] [i - 2sin 2 (Mcos( 2M 


48v/2tt 3 


2 2 2 
i 2 m 7 , t m 7 r TT 

log —f-41og—f+5- — 
mi mi 3 


(7) 


As in the calculation of the W —> pvcj) decay, we have 
neglected the muon mass, and we have expanded the Z’s 
decay width in Q to leading order in m^/mz- These 
steps are motivated by the arguments given in the para¬ 
graph following (|5|. 


III. SCALAR THEORY 

We also consider a scalar theory which is well behaved 
without the addition of any shadow particles. The inter¬ 
action Lagrangian is 


£mt,S = 4>S 
+ <t>P 


Cstpfifpu. + C'gV’pV’p 


( 8 ) 


+ CPpi/jp^ipp 


where <ps is the scalar field, 4 >p is the pseudo-scalar field 
where TO 0 s = , and the (7’s (with corresponding su¬ 

perscripts and subscripts) are the corresponding coupling 
strengths. In this section it is understood that 4> refers 
to either <ps or </>p . 

As with the vector theory, we again consider the con¬ 
straint due to the branching ratio of W —> pvips plus 
W -A g,v(f)p. The decay amplitude for both scalar and 
psuedoscalar cases is given by the Feynman diagram in 
Fig. |5l 
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FIG. 5: W -> nv</> 


The result for the decay width is 


scalar 


G F m 3 w [(C%) 2 + (Ctf 
96-\/27r 3 




m<f, (MeV) 


where, as with the vector case, we have taken the muon 
and neutrino mass to be zero. 

As with the vector theory, this logarithmic dependence 
in opens the parameter range in Cg necessary to solve 
the proton radius problem. This is seen in the constraint 
plot Fig. [6] 

IV. CONCLUSION 

We have considered some of the constraints faced by 
proton radius puzzle explanations that hypothesize a new 
muon-specific force. 

In particular, we have considered some of the necessary 
features required for a vector mediator to be embedded 
into a larger renormalizable model. Without the em¬ 
bedding, amplitudes grow quickly with energy, and the 
energy in W decay is enough to cause trouble. The ra¬ 
diative corrections to W —> /jlv would be larger than the 
uncertainties in the measured and well calculated W de¬ 
cay rate, using couplings that could explain the proton 
radius puzzle. However, with a theory which also con¬ 
tains Yang-Mills couplings, the unphysical energy depen¬ 
dence is ameliorated and the W width constraint causes 
no problem for BSM models solving the proton radius 
puzzle. 

The tightest remaining constraint follows from some 
measured energy splittings in muonic Mg and Si, which 
are well calculated in the standard model, and new 
physics cannot cause a deviation larger than the exper¬ 
imental uncertainty. Though this constraint is tight, it 
does eliminate the muon-specific interaction at the level 
required for the proton radius problem. 

Thus, it is fair to say that a theory involving a new 
vector (and axial vector) boson, <fi, along with vector 
(and axial vector) couplings to the muon and proton, 
as described by the Lagrangian in Eq. ([I]), gives a viable 
representation of the results for some to-be-determined 
complete solution to the proton radius puzzle. Such a 
theory necessitates the existence of an additional W-like 
boson, the shadow W denoted W a , which fixes the other¬ 


FIG. 6: The parameter space necessary to satisfy experimen¬ 
tal constraints on the scalar coupling Cg. The green band, 
outlined by solid lines, is the constraint on Cg necessary to 
solve the proton radius problem (±2cr). The shaded red re¬ 
gion is the restricted region on Cg due to the requirement 
that the branching ratio for W goes to g,v(j>s + g-v<f)p must be 
less than 4% under the assumption that Cp solves the muonic 
g — 2 problem. The shaded orange region is the restricted re¬ 
gion on Cg due to energy splittings in muonic Mg and Si at 
2 <j. The dashed green line is the constraint on Cp necessary 
to solve both the proton radius and muonic g — 2 discrepancies 
(for clarity, only the central value is shown). 

wise poor high energy behavior in radiative corrections to 
W —> decay first noticed by Karshenboim et al. [TS] . 

The arguable success of the vector and axial vector 
theory, or the scalar and pseudoscalar analog, motivates 
finding a suitable complete embedding into the Standard 
Model. Such an embedding must be renormablizable. 
The addition of a shadow W and corresponding triple 
boson interactions is one step in this direction (as pointed 
out by Llewellyn Smith ED). 

Finding such an embedding may necessitate new in¬ 
teractions which may lead to additional constraints. For 
example, one must consider the process by which the 
new bosons acquire their mass. We also expect loop cor¬ 
rections to the self energy of the W. Corrections with 
loops of new particles correlate with contributions to 
the Peskin-Takeuchi parameters [25H25| . Though these 
parameters are only defined under certain conditions, 
e.g., for BSM theories in which new particles are heavy, 
one can find analogous parameters for theories in which, 
as in our case, the new boson is light. Contributions to 
these parameters will be calculated and compared against 
experimental constraints. 

It is also relevant to note that the addition of a new 
muonic interaction necessitates new radiative corrections 
beyond just W decay. For example, one may also con¬ 
sider muon decay, pion decay, and the like. Such pro¬ 
cesses are at a significantly lighter mass scale compared 










6 


to that of W decay and thus are less likely to cause trou¬ 
ble. We plan to asses if this is true in the near future. 
Furthermore, changes to the rate of muon pair produc¬ 
tion from proton collider experiments may provide fur¬ 
ther constraints on the theory. 

These are just some examples of how a vector boson 
theory may be more restricted than Fig.[3]appears. How¬ 
ever, successes so far point to the still possible viability 
of these theories as solutions to the proton radius (and 
muonic g — 2) puzzle. Thus, continuation of research 
down this track is well motivated. 
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